Three-body breakup within the fully discretized Faddeev equations 
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A novel approach is developed to find the three-body breakup amplitudes and cross sections within 
the modified Faddeev equation framework. The method is based on the lattice-like discretization of 
the three-body continuum with a three-body stationary wave-packet basis in momentum space. The 
approach makes it possible to simplify drastically all the three- and few-body breakup calculations 
due to discrete wave-packet representations for the few-body continuum and simultaneous lattice 
representation for all the scattering operators entering the integral equation kernels. As a result, the 
few-body breakup can be treated as a particular case of multi-channel scattering in which part of the 
channels represents the true few-body continuum states. As an illustration for the novel approach, 
an accurate calculations for the three-body breakup process n+d —¥ n+n+p with non-local and local 
NN interactions are calculated. The results obtained reproduce nicely the benchmark calculation 
results using the traditional Faddeev scheme which requires much more tedious and time-consuming 
calculations. 
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I. INTRODUCTION 



The last decades have inaugurated great success in 
precise ab-initio calculations for few-body scattering pro- 
cesses [J-01- These calculations made it possible to de- 
scribe accurately the results of numerous recent experi- 
ments on elastic nd scattering at energies up to 350 MeV 
and also the three-body breakup n + d — > n + n+p at low 
and moderate energies E n ~ 10-^30 MeV. However some 
problems remain unsettled even at such low energies. 
These are the so-called A y -puzz\e (as well as other puzzles 
for various tensor and vector analyzing powers) in elastic 
scattering, the problems with an adequate description 
of the pairwise 1 5'o-channel contribution to three-body 
breakup at low energies Q and breakup cross section in 
some particular three-particle configurations such as the 
quasi-free scattering Q and the space-star [l(| configura- 
tions. The most plausible reason for the visible discrep- 
ancies with experimental data in this area is likely not 
insufficient accuracy of numerical calculations but rather 
some deficiency in the input 2N- and 3iV-interactions. 
At the same time, the progress in the field of precise few- 
nucleon calculations, particularly in testing of new mod- 
els for 37V-interactions, is restrained strongly by a high 
number of very complicated few-nucleon calculations, es- 
pecially above the three-body threshold. Because of these 
complications of traditional computational schemes for 
the direct solution of the Faddeev- Yakubovsky equa- 



*Electronic address: rubtsova-olga@yandex.ru 

^Electronic address: pomeran@nucl-th.sinp.msu 

* Electronic address: kukulin@nucl-th.sinp.msu 

§ Electronic address: faessler@uni-tuebingen.de 



u.ru| 
.ru| 



tions, one observes a rise of interest in recent years in 
alternative approaches [IlT - fl3] ] to calculate the scatter- 
ing observables by simpler methods. 

Among such alternative approaches one can note a 
preference for the so called L 2 -methods. These meth- 
ods are based on expansions of the scattering solution 
into a basis of square- integrable functions [1J-[21|. Such 
L2-methods proved to be very well suited and quite 
efficient for numerous applications. One of the most 
successful approaches of this type is the Continuum- 
Discretized Coupled- Channel (CDCC) method developed 
three decades ago for treatment of breakup processes in 
direct nuclear reactions fl9U22j ]. The CDCC approach 
in its traditional form was unable to treat other channels 
than elastic scattering and projectile (or target) breakup. 
Recently a few groups generalized the traditional CDCC 
approach to scattering of three-fragment projectiles by 
a stable target [2l|. However this generalized approach 
can be considered as a hybrid method: L2 discretization 
of inner motion in the three-body projectile and the tra- 
ditional treatment of a coupled-channel problem 

On the other hand, the present authors have developed 
some alternative L2-technique f23rf26j which is based on 
the idea of complete continuum discretization with a spe- 
cial stationary wave-packet basis in momentum space 
(three-body lattice basis). The basic distinction of such 
an approach from the traditional CDCC-scheme for the 
three-body systems is that the wave-packet approach is 
dealing with a full discretization of the three-body con- 
tinuum. In other words, the discretization on both Jacobi 
coordinates is used here rather than the discretization on 
the alone coordinate of the projectile inner motion as in 
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the CDCC approach 1 . 

Our approach with the global discretization over all 
valence coordinates leads immediately to a few impor- 
tant advantages in the accurate treatment of few-body 
scattering. Among those the following are the most im- 
portant: 

(i) The few-body scattering problem is consistently 
formulated in a Hilbert space of three-body nor- 
malized states, similarly to the bound-state prob- 
lem. 

(ii) The approach employs the integral equation frame- 
work of scattering theory instead of the differential 
equation approach (e.g. in the CDCC) where the 
boundary conditions in few-body scattering chan- 
nels are not easy to formulate, especially in terms 
of the L2 basis used. Contrary to this, the integral 
equation formulation allows to avoid any explicit 
account of the boundary conditions. 

(iii) When working within the wave-packet formula- 
tion of the scattering problem one can derive ex- 
plicit formulas for some scattering operators (e.g. 
channel resolvents). Such fully analytical finite- 
dimensional approximations, being substituted into 
integral equation kernels, lead immediately to their 
algebraic matrix analogues. Thus, our final equa- 
tions are simple matrix linear equations with regu- 
lar matrices. 

In our previous works [13, HH we have demonstrated 
how to find the elastic 2 — » 2 scattering amplitudes in lat- 
tice representation. In the present paper, we generalize 
the technique to the three-body breakup treatment. So 
we present here the complete formalism for determina- 
tion of three-body breakup amplitudes. It is important 
to emphasize in this connection that the accurate treat- 
ment of three-body breakup within the Li type approach 
is much less obvious than that of elastic ones and thus 
requires some additional delicate theoretical studies. In 
particular, the matrix elements in the breakup amplitude 
are not truncated over all spatial coordinates (in contrast 
to the elastic and rearrangement amplitudes), so the va- 
lidity of the L2 scheme in the treatment of the breakup 
processes should be studied carefully. As some substanti- 
ation for such approach one can consider the three-body 
breakup calculations within the CDCC-approach where 
the discretization of the continuum in the projectile in- 
ner subHamiltonian has been used for the description of 



the breakup amplitudes [22|. So, the natural general- 
ization of such a partial continuum discretization to the 
case of full three- and few-body continuum within the 
Faddeev equation approach is an important next step. 
Moreover, this fully discretized approach studied in the 
present work allows to simplify drastically all calculations 
and makes it more universal and elegant. 

The present work has the following structure. In the 
section II, a three-body lattice-like free wave-packet ba- 
sis is described in detail together with a similar basis for 
the channel Hamiltonian. Here we also discuss the prop- 
erties of these bases. The complete formalism for elastic 
scattering and breakup, as applied to the nd-system in 
the packet representation is presented in Section III. In 
the Section IV, a few useful numerical illustrations and 
their comparison with the standard Faddeev benchmark 
calculations are given. Our results are summarized in the 
Section V. For the sake of convenience for the reader we 
add three appendices. In Appendix A we describe the 
detailed scheme for calculation of the three-body over- 
lap matrix in the three-body lattice basis for recoupling 
of different Jacobi coordinates. In Appendix B we give 
the convenient wave-packet formalism for the solution of 
three-body scattering problem with separable pairwise 
interactions. In the last Appendix C we discuss some 
features of our numerical calculations. 



II. LATTICE REPRESENTATION FOR THE 
THREE-BODY CONTINUUM 

We consider here the problem of a scattering of three 
identical particles 1, 2 and 3 (nucleons) with mass m, 
interacting via pairwise short-range potentials v a {a — 
1,2,3). It is convenient to use three Jacobi momentum 
sets (pa,q a ) corresponding to three channel Hamiltoni- 
ans H a (a — 1,2,3) which define the asymptotic states 
of the system. For example, the channel Hamiltonian Hi 
has the form of the direct sum of two-body subHamilto- 
nians 

H x = h x ® hi (1) 

where subHamiltonian hi describes NN subsystem con- 
sisting of particles 2 and 3 with interaction vi and sub- 
Hamiltonian h 01 corresponds to the free motion of nu- 
cleon 1 relative to center of mass of the subsystem {23}. 

As we study the identical particle system, we will omit 
further, where it is possible, the Jacobi index a. 



It is interesting to add that similar idea of global three-body 
discretization in a momentum space has been proposed earlier 
|l. r | within the pseudostate extension of the coupled-reaction- 
channel method. The author solved as an illustration of the 
approach the simple model problem of 2— >2 scattering and also 
the breakup 2 — ¥ 3 process using the Laguerre polynomial basis. 
Unfortunately this prospective approach has not received any 
further development. 



A. The two-body free wave-packet states 

We start from the free-motion three-body Hamiltonian 
defined in the given Jacobi momentum set (p,q) 

H = h ® hi (2) 
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where the subHamiltonian ho defines the free motion of 
two nucleons with the relative momentum p and the sub- 
Hamiltonian Kq defines the free motion of the third nu- 
cleon with the momentum q relative to the center of mass 
of given N N subsystem. 

Now we will construct our three-body L2 basis using 
discretization of the continua of the two above subHamil- 
tonians. In doing this we will employ the complete sets 
of continuum wave functions \p) and \q) (for every partial 
wave) normalized according to the conditions: 

{p\p')=5(p-p'), (q\q')=S(q-q'). (3) 

When discretizing, we truncate the continuum of ho 
and by maximal values e ma x and £ m ax respectively, 
so that the continuous spectra above these values can 
be neglected. Further, the selected energy regions 
[0,e m ax] and [Ojfmax] are divided onto non-overlapping 
bins {{ei-i^iltU} and {{ S j-U S j\f=i}- Sucn energy bins 
correspond to momentum bins [pi-i,pi] and [qj—%,qj], so 
that the end-points of both sets are interrelated by con- 
ventional formulas pi = ^/mel and qj = ^m£j. To 
further simplify the notation, we will denote the inter- 
vals in the variable p (both the energy and momentum 
ones) as ®i and in the variable q as We use also the 
following notations for the widths of the corresponding 
momentum intervals: 

di =Pi- Pi-i, dj = q 3 - qj-i (4) 

Now let's define a set of free stationary wave packets 
(WPs) as integrals of the plane waves (corresponding to 
the free motion) over the above momentum bins for both 
subHamiltonians 2 : 

\pi) = -L= [ dpf(p)\ P ), i = l,..., M, (5) 

|<b-> = -7== / dqw{q)\q), j = l,...,N. (6) 

where f(p) and w(q) are some known weight functions 
and Ai and Bj are normalization factors, directly related 
to the weight functions 

Ai = [ dp\f{p)\\ Bj = f dq\w{q)\\ (7) 

so that the WP states are normalized to unity: 

(pi\pi') =S U ', ((\Mj')=S jf . (8) 

It is important to stress, these WP states belong to a 
Hilbert space (similarly to the bound state functions) 
and WP functions are square-integrable: in configura- 
tion space they vanish at infinity in contrast to the initial 



plane waves. But in the relevant restricted range of con- 
figuration space the WP states still resemble quite closely 
the exact scattering states taken at the bin center energy 
(or momentum) |24j |. 

The sets of such WP states \pi)f£i and \qj)f = i form an 
orthonormalized bases in Hilbert space which can be used 
as normal L2 bases, e.g. also for variational calculations. 

In our previous papers [2J, |25| we have discussed the 
properties of WP's in detail. A distinctive feature of WP 
bases is that the matrices of the subHamiltonians found 
in such bases are diagonal: 

(pi\ho\pi>) = e*S u >, (qMqj') = £*8 jr , (9) 

where values e* and E* arc defined via corresponding end- 
points of bins Di and S)j [24[. The most useful property 
of WPs is that the matrices of the resolvents go(e) — 
[e + iO — ho]^ 1 and g\{£) = [£ + iO — ^o] -1 are diagonal 
in the corresponding WP bases and their elements have 
explicit analytical forms [24| . 

Different choices of weight functions lead to different 
sets of WPs. In practical calculations in this work we 
use the momentum wave packets with the unit weight 
functions: 

f(p) = 1, At = d h w(q) = 1, Bj = dj. (10) 

It is easy to see that the overlap of such free momentum 
WP with a plane wave, i.e. the momentum representa- 
tion of packet state itself has the form: 

<Mh> = (ID 

where we have introduced a function d(p G 35 j) which is 
equal to unity if the momentum p belongs to the interval 
[pi~i,Pi] and vanishes in the other case. So, the wave 
packet |pi) takes a form of simple step-like function in 
the momentum representation. 

The sets of the constructed free WP states can be used 
to find two-body bound-states and also to solve a two- 
body scattering problem, e.g. for finding the two-body 
off-shell i-matrix [2J|. In the present paper we will use 
these two-body L 2 bases to construct three-body WPs 
for solution of the three-body scattering problem. 

B. Three-body lattice basis and the permutation 
matrix 

Three-body wave packet states are built as direct prod- 
ucts of two-body ones. However, here one should take 
into account the spin and angular parts of the functions. 
The total three-body WP basis function can be written 
as: 

\XW > ) = \p?)®\q*)\ a ,l3:r), (12) 



2 Below we will use the Gothic letters to denote objects (wave 
functions and operators) in the WP subspace. 



where \a) is the spin-angular state of the NN pair, 
is the spin-angular state of third nucleon, while |T) is 
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the set of the three-body quantum numbers. The state 
(fT2")l is a WP analog of the exact state of the three-body 
continuum \p, q)\a, j3 : T) for the free Hamiltonian Ho- 

The properties of such three-body WP's are very simi- 
lar to those of two-body wave packets 24J . In particular, 
the matrix of the three-body free Hamiltonian Ho and 
its resolvent Gq(E) = [E + iO — Hq]^ 1 are diagonal in 
the so constructed basis. In other words, such a WP ba- 
sis defines an "eigen" wave-packet subspace for the free 
three-body Hamiltonian Ho. 

Since the basis functions are the products of both 
step-like functions in variables p and q, the solution of 
the three-body scattering problem in such a basis corre- 
sponds to a formulation of the scattering problem on a 
two-dimensional momentum lattice. Therefore we will re- 
fer to such a basis as a lattice basis. Let us denote the two- 
dimensional bins (i.e. the lattice cells) as 55 jj = 55i<8S5j. 
In the few-body case, the lattice basis functions are con- 
structed as direct products of the two-body free WPs, so 
the basis space corresponds to a multi-dimensional lat- 
tice. 

In principle, using the above lattice basis, one can solve 
a general three-body scattering problem by projecting all 
the scattering operators onto such a basis. In particular, 
the matrix of the three-body free resolvent Go can be 
expressed in the above lattice representation fully ana- 
lytically Hf|. 

Let us consider the particle permutation operator P 
which enters in the Faddeev equation for three identical 
particles and is defined as 



P — P12P23 + PviPiz- 



(13) 



The matrix of the operator P in the lattice basis cor- 
responds to the overlap between basis functions defined 
in different Jacobi sets: 



(X^\P\Xf,f) = (X^(1)\X^ (2)), (14) 

where the argument 1 (or 2) in the basis functions means 
a corresponding Jacobi set. Such matrix element can be 
calculated with the definition of the basis functions in 
momentum space (|11|) : 



Tap' 



of the permutation operator matrix in the WP basis are 
finite. 

Using the above "packetting" procedure and the hy- 
perspherical momentum coordinates, the calculation of 
the matrix element in Eq. (|15[) can be done using only a 
one-dimensional numerical integration over the hypermo- 

mentum K = ^Jp 2 + |g 2 . The technique of this calcula- 
tion for the s-wave basis functions is given in Appendix 
A of the present paper. The generalization for higher 
partial waves is straightforward. 

It should be emphasized here that the fixed lattice-like 
form for the permutation operator matrix makes it possi- 
ble to avoid the complicated and time consuming multi- 
dimensional interpolations of the current solution when 
solving the Faddeev equations (in momentum space) by 
iterations in conventional approach [l|, Q . Such numerous 
multi-dimensional interpolations at each iteration step 
take a big portion of the computational time in prac- 
tical numerical procedure. When solving the four-body 
Yakubovsky equations the dimension for these interpo- 
lations increase and thus the computational efforts get 
even higher. So, avoiding the very numerous multidimen- 
sional interpolations in each step of the iterations leads to 
a tremendous acceleration for all three-body calculations 
in momentum space. 

Thus, the two-dimensional momentum lattice basis 
constructed above can be applied directly to solving the 
Faddeev equations for the conventional transition oper- 
ator U. However, by using the very convenient form 
for the spectral representation of the resolvent operators 
in the WP basis one can employ some alternative (but 
equivalent) form of the Faddeev equation, which makes 
it possible to avoid the time-consuming calculation of 
the fully off-shell t-matrix at many energies (which re- 
quires to solve very often the Lippmann-Schwinger equa- 
tions at every energy and for different spin-orbit chan- 
nels) and to replace it by calculating the resolvent of the 
NN subHamiltonian hi in the corresponding scattering 
WP-representation. The latter can be made easily by 
straightforward one-fold diagonalization of the hi sub- 
Hamiltonian matrix. 



l0 ]^' - {x^\p\xlf) = 



dpdq / dp dq x 

Si 3 ., 



, (15) 



di d> j ' 



where the prime at the lattice cell 35^-, indicates that 
the cell belongs to the other Jacobi set while the 
P£p a 'p'(p,q,p'q') is the kernel of particle permutation 
operator in a momentum space. This kernel, as is well 
known [27j . is proportional to the product of the Dirac 
delta and Heaviside theta functions. However, due to 
"packetting" (i.e integration over momentum bins in 
Eq. (fPo]) ) these singularities get averaged over the cells 
of the momentum lattice and, as a result, the elements 



C. The scattering WPs for the subHamiltonian hi 



As has been demonstrated earlier (24| - |26| the station- 
ary wave packets can be built not only for free Hamilto- 
nians but also for perturbed two-body hi = ho + vi and 
three-body channel Hamiltonians Hi . 

In case of the hi subHamiltonian, its continua [0, e max ] 
for every spin-angular configuration a are divided into 
separate bins e k\h=ii an d one can build the scat- 

tering wave packet for every such bin = in 
the form 




(16) 
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i.e. as an integral over the exact scattering wave function 
\4>p) on the energy interval A£. Here we use the unit 
weight function and D% is the width of interval Ag . 

It is easy to show that such packet states have the same 
properties with respect to their "eigen" Hamiltonian hi 
as free WPs with respect to the free Hamiltonian ho. The 
only difference is that the set of scattering WP's should 
be accomplished with the bound state functions of h\ (if 
they exist). Jointly with the possible bound state wave 
functions for the h\ subHamiltonian, the scattering WP's 
form an orthonormalized basis in which both the matrix 
of the Hamiltonian h\ and the matrix of its resolvent 
gi(e) = [e + iO — /ii] -1 are diagonal (2^ |. 

The projection properties for the WP of hi will be 
similar to those for ho viz.: 

U a, a , gfr € Ag) 

\<Pp\ z k) = rj^ ■ ( 17 ) 

V k 

D. Pseudostates as approximations for scattering 
WPs 

At first glance it may appear that the exact scattering 
WP basis is useless because its construction would re- 
quire the knowledge of exact scattering wave functions of 
the Hamiltonian h\. However, as has been demonstrated 
[24j . the properties of the exact scattering WP's for h\ 
are quite similar to those of respective pseudostates ob- 
tained by the diagonalization of the Hamiltonian matrix 
in some complete L2-basis. So that, in actual calculations 
one can replace the set of WP's \z%) by the set of respec- 
tive pseudostates [24[ . Such an L2-basis can be used as 
a very good approximation for the free WP-basis ([5]). As 
a result of such Hamiltonian matrix diagonalization, one 
gets a set of pseudostates 

M 

l^> = £0&ltf>, k = l,...,M, (18) 

2=1 

together with a set of their eigenvalues e%* . 

In this paper, we restrict ourselves to s-wave spin- 
dependent pair interactions only. We assume that there 
is a single bound state \z ) (deuteron) with binding en- 
ergy £g in the AW spin-triplet channel and there are no 
bound states in the A A spin-singlet channel. 

In case of s-wave scattering with s-wave A A interac- 
tions, the indices a, (3 and T in Eq. (fl~2j) include only the 
spin quantum numbers. So, below we will use the value 
of the spin of AA-pair s = 0, 1 instead of index a, while 
the index /3, which indicates the spin value of the third 
nucleon (i.e. 5), will be omitted everywhere. The index 
r defining the set of quantum numbers for three-body 
states is reduced to the total spin of the three-body sys- 
tem S = I, I which in the s-wave case is equal to the 
total angular momentum of the system and therefore is 
conserved. 

After the above diagonalization in the spin-triplet 
channel (s = 1) one gets a set of pseudostate functions, 



the first of which \z\) w \zq) with the energy £q, is an 
approximation for the deuteron wave function, while the 
other M — 1 pseudostates with energies e\* are localized 
in the continuum spectrum and correspond to scattering 
WP states for hi. In the spin-singlet channel there are no 
A A bound states, so that all functions \zf*) in Eq. (fT8|) 
are approximated by scattering WP's. It is important to 
note that as a by-product of our diagonalization proce- 
dure one gets simultaneously the discrete representation 
for A A partial phase shifts S s (el*) for all pseudo-states 
energies (i.e. in one step!) - see the detail in Ref. p8|. 

Since the free WP basis functions (in the momen- 
tum space) are step-like functions, the momentum de- 
pendence of all functions expressed via such a basis have 
a histogram-like form. An example of the momentum 
dependence for the bound state (deuteron) function in 
such step-like basis in comparison with the exact func- 
tion for Yamaguchi triplet AA potential (see Appendix 
B) is displayed in Fig. [TJ 
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FIG. 1: (Color online) Comparison of the exact deuteron wave 
function obtained in the momentum space for the Yamaguchi 
potential (dashed curve) with its approximation in the lattice 
basis (solid line). 

The Fig. [2] displays the functions of two pseudostates 
(with k = 4, 8) obtained in the lattice basis in compari- 
son with the corresponding exact scattering wave packets 
which can be calculated exactly for the separable Yam- 
aguchi potential. It is interesting to see that although 
functions of the exact scattering WP's (|TB)) (dashed lines 
in the figure) have the logarithmic singularities at the 
boundaries of the "on-shell" interval (i.e. the one which 
the state energy belongs to) they are square-integrable 
as well as the free-motion WP's. 

It is clear from the comparison that the pseudostates 
composed from step-like wave packets reproduce very 
well the structure of the exact scattering wave packets 
"on average" . 

Now having at our disposal the two-body bases for 
subHamiltonians hi and h$ one can construct the three- 
body WP basis for the channel Hamiltonian Hi which 
defines the asymptotic motion in the system. 
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FIG. 2: (Color online) The functions of pseudostates (k — 
4, 8) obtained in the lattice basis (solid lines) in comparison 
with exact scattering packets (dashed lines) for the AW spin- 
triplet Yamaguchi potential. 



E. Construction of three-body WP basis for the 
channel Hamiltonian 

The three-body WP states corresponding to the chan- 
nel Hamiltonian H\ can be defined similar as the WP- 
states for the three-body free Hamiltonian Ho, i.e. as di- 
rect products of two-body WP states for Hq and hi sub- 
Hamiltonians (jointly with the bound-state) multiplied 
by the spin functions of the system: 

\Zf») = \4)®\ qj )\s±:Z), (19) 

where s and £ are the NN subsystem and the total three- 
body spins correspondingly 3 . 

When using the above pseudostate approximation, 
these three-body states, similarly to two-body scattering 
WPs, are related to the three-body lattice basis states by 
a simple rotation transformation (similar to Eq. (|18l) ): 

M 

\zfn=J2°ki\ x ^ s )- (2°) 

i=l 

Hence, starting from the free WP bases for every pair 
subsystem one gets a set of basis states both for the three- 
body free Hq and the channel H\ Hamiltonians. The 
basis defined in Eq. (1191) defines an "eigen" WP-subspace 
for the channel Hamiltonian H\. 

This allows us to construct an analytical finite- 
dimensional approximation for the channel resolvent 



Gi(E) = [E + i0 - Hi]' 1 . Indeed, the exact three-body 
channel resolvent is the convolution of the two-body sub- 
resolvents gi (e) and g$(£)- 

i r°° 

Gi(i?) = — y deg 1 (e)g 1 (E-e). (21) 

Using further the spectral expansions for the two-body 
resolvents and integrating over s, one gets an explicit ex- 
pression for the channel resolvent Gi as a sum of two 
terms G X {E) = Gf c (E) + G^ C {E). Here the bound- 
continuum part Gf c (E) is the spectral sum over the 
three-body states corresponding to the free motion of 
the deuteron relatively to the third nucleon. So, the 
imaginary part of Gf c (E) is related to a discontinuity 
on the two-body cut of the Riemann surface of the three- 
body energy E. The continuum-continuum part Gf C (E) 
of the channel resolvent includes the channel three-body 
states with the iViV pair with interaction in the contin- 
uum and ImGf c (E) is defined by a discontinuity on the 
three-body cut on the energy surface (see the details in 
Ref. 0). 

Projecting the exact channel resolvent onto the three- 
body channel WP basis defined in Eq. dl~9l , one can find 
analytical formulas for the matrix elements of the Gi 
operator. The respective matrix is diagonal in all wave- 
packet and spin indices: 

(ZftldWlzft,) = 5 kk ,S 3f 6 ss ,Gf;(E). (22) 

Here the diagonal matrix elements G^j 5 (E) are defined 
as integrals over the respective momentum bins and de- 
pend in general on the spectrum partition parameters 
(i.e. the pi and qj values) and the total energy E only. 
They do not depend explicitly on the interaction poten- 
tial vi . If the solution of scattering equations in the finite- 
dimensional WP basis converges with increasing the ba- 
sis dimension, the final result turns out to be independent 
on the particular spectral partition parameters. We have 
found [261 ] the explicit formulas for the resolvent matrix 
elements (|2"2"|) when one uses the energy WP's 4 , i.e. WP's 
with the weight functions f(p) = */p, w(q) = ^fq. 

The representation (1221) for the channel resolvent is the 
basic expression for our wave-packet approach, since it 
gives explicit analytical formulas for the three-body resol- 
vent and thus it allows to simplify drastically the solution 
of general three-body scattering problem. This expres- 
sion can be used directly to solve the finite-dimensional 
analog of the Faddeev equations for the three compo- 
nents of the total scattering wave function ,25]. Al- 
ternatively the very convenient representation (|22j) can 



3 We consider here the three-body states with total isospin T = i 
only. Since in case of s-wave NN pairwise interaction the spin- 
and isospin quantum numbers are interrelated uniquely by the 
Pauli principle we can omit further the isospin parts of the wave 
functions and corresponding quantum numbers. 



4 The matrix elements of the three-body channel resolvent take a 
simple analytical form in the WP basis constructed from the con- 
tinuum wave functions normalized to the delta-function on the 
energy (the energy WPs). For finding the resolvent matrix ele- 
ments with WP's with various weight functions one uses renor- 
malization factors for a transition from the given wave packet 
states to the energy ones. 
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also be used to solve some particular three-body scatter- 
ing problems using the three-body Lippmann-Schwinger 
equations 24]. 



III. SOLUTION OF nd SCATTERING PROBLEM 

Now let us proceed with solving the nd elastic and 
breakup scattering problems. 



A. The elastic and breakup rid scattering 
amplitudes 

The elastic scattering observables can be found from 
the Faddeev equation (FE) for the transition operator U, 
e.g. in the form 



U = PGn 



PtG U, 



(23) 



where t is two-body 7j-matrix in three-body space and 
P is the particle permutation operator. The equivalent 
form of FE for the transition operator U has the form: 



U = Pv 1 + Pv 1 G 1 U, 



(24) 



where G\ is the resolvent of the channel Hamiltonian 
H\. Since tGo = v\G\ the operators U and U coincide 
on-shell and half-shell. 

Since the determination of the off-shell channel resol- 
vent in three-body space is a rather time consuming so- 
lution of the FE in the form (|24|) . it is very seldom em- 
ployed for practical solutions. Actually a similar form of 
the equations is associated with formalisms of the con- 
figuration space Faddeev equations, where quite differ- 
ent numerical approaches have been employed [29], |3(| 
than for the momentum space FE. However, since in the 
lattice approach one has explicit analytical formulas for 
the three-body channel resolvent G\ the form (pH)) of FE 
turns out to be very appropriate for the numerical solu- 
tion in a WP basis. 

The elastic nd scattering amplitude (for a given value 
of total spin S) can be defined as matrix element of 
the solution of the Eq. (|24l) taken in the initial state 
1^0)30, E): 



A* 



fob) 



2 m 

3 go 



z ,qo,^\U\z ,q ,T,} 



(25) 



The breakup amplitude for one Faddeev component 
of the three-body wave function (the so-called single- 
component amplitude) can be found from the elastic 
transition operator U after applying the operator tGo 
from the left: 



T^(p,q) 



(p L g 1 E£|tgot7]^o 1 go 1 E) 
pqqo 



(26) 



To obtain the differential breakup cross sections, the total 
breakup amplitudes can be found by the contributions of 
all three single-component amplitudes. 



Thus, we change the conventional treatment of the 
breakup process [l[ and consider the three-body asymp- 
totic states as scattering states for the channel Hamilto- 
nian Hi rather than as the states of the three-body con- 
tinuum for the free Hamiltonian Hq . It is a quite natural 
when treating the elastic scattering amplitudes because 
the initial state wave functions correspond to the chan- 
nel Hamiltonian Hi. In full analogy with this, one can 
treat the deuteron breakup as its excitation into a con- 
tinuum A.ZV-state in the two-body subsystem governed 
by the hi subHamiltonian 5 . As was already indicated 
above such a treatment of breakup processes is close to 
the configuration space approach. 

In fact, when solving the three-body Faddeev equa- 
tions in the configuration space [29l - [3lj one finds the 
breakup amplitude A(6) which determines the asymp- 
totic behavior of the three-body wave function in hy- 

perspherical coordinates p = 

arctan(-7=-) as follows 



^y 2 and ■& — 



^ y) ^WpW-' 



K = 



P 



-q 

A 1 



(27) 



where x and y are two Jacobi coordinates and K is the 
hypermomentum. 

This breakup amplitude is defined for every spin- 
angular configuration and interrelated to the partial 
single-component breakup amplitudes (|26p by the follow- 
ing formula 

3v3 2p 

(28) 

where 9 is the hyperangle in momentum space. 

Now if one transforms the formulas for the breakup 
amplitudes A from [2{| to the integral form one receives 
the following definition for the breakup amplitudes in the 
momentum hyperspherical representation 

4-7777? K 4 

A^ s (9) = — —eW^^vji^W^ (29 ) 



3^ pq 



where \4>p ) is scattering function for the Hamiltonian 
hi corresponding to the outgoing boundary condition. 
These functions are distinguished from the real-valued 
functions \4>p) used in our approach by only a phase fac- 
tor: 



(30) 



where S s (p) is the s-wave phase shift of the NN scattering 
in the channel with spin s. 



5 It is of interest to remark that while just such a scheme has been 
used in the CDCC treatment of breakup processes |22| , in the 
Faddeev approach the final states used for the breakup treatment 
are the free three-body states. 
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Using formulas ((29)) and (|28|) one can derive an alter- 
native to formula (|26[) for the single- component breakup 
amplitude via the scattering functions of the channel 
Hamiltonian Hi 



pqqo 



(31) 



Summarizing this derivation one can conclude that the 
breakup amplitudes can be defined quite similar to a ma- 
trix element for the elastic scattering transition operator 
U with replacement of the the TVTV bound-state wave- 
function with the exact scattering functions for the TV TV 
subHamiltonian. 

Having now the required representations for both the 
elastic and breakup amplitudes, we will proceed in solv- 
ing the Faddeev equation in "eigen" WP subspace of the 
channel Hamiltonian Hi . 



Solution of the Faddeev equation in the 
three-body WP basis 



(1351) can be found analytically for a wide variety of the 
potential forms. 

An important ingredient of our lattice approach pre- 
sented here is the representation of the permutation op- 
erator P as an overlap matrix P between the channel 
WP basis functions for different sets of the Jacobi co- 
ordinates. Using further the approximation (fT5|) for the 
scattering wave packets \z k ), these matrix elements can 
be expressed through the overlap matrix P° for the free 
lattice basis function of Eq. (|T5|) with the help of the 
rotation matrices O 



(Z%\P\Z%>) 



(36) 



Now let's replace the exact operator U in the formula 
for the elastic nd scattering amplitude (|25|l with its lat- 
tice counterpart and employ further the projection 
rule for the free WP states (fTTj). Then one gets that the 
on-shell elastic amplitude in the wave-packet representa- 
tion can be calculated as a a diagonal (on-shell) matrix 
element of the U-matrix: 



In our wave-packet approach, all the operators in 
Eq. I|24p are projected onto a three-body wave-packet 
basis corresponding to the channel Hamiltonian H±. In 
other word, every operator, e.g. U, is replaced with its 
finite-dimensional WP representation: 



il S = EE \Z%)(Z%\U\Z$j,)(Z%,\. 



(32) 



s,kj s' ,k' j' 



Finally one gets the matrix analog for the Eq. (|24|) (for 
the given value of S) 



U = PVi + PViGiU. 



(33) 



Here Vi and Gi are the matrices of the pair interaction 
and the channel resolvent respectively, the matrix ele- 
ments of which can be found in an explicit form. 

The matrix Vi of the potential v\ is diagonal in the 
indices j, j' of the wave-packet basis ((HJ) for the free sub- 
Hamiltonian h° and has the block form: 



]kik>j> = 8ij'8.A4\v'iW)- 



(34) 



These matrix elements do not depend on in index j and 
can be written with the usage of the rotation matrix O 
defined in Eq. (JTHJ) as: 



In the last expression, the potential matrix elements in 
the free WP basis are used which have the form 



(pJKIp*') 



i 



y/didi 



dp / dp' vl(p,p') (35) 



where v{(p,p') is the momentum representation for the 
interaction potential. It implies that the matrix elements 



2m 



3<?o 



E = 



'Jo 



1 3 

2' 2' 



(37) 



where l-^ojo) is the WP basis state corresponding to the 
initial state: index denotes the bound state of the TV TV 
pair (deutcron) and index j denotes the "on-shell" q-bin 

i)j with the on-shell momentum qo = J ^m(E — £q): 

9o € If jo. 

It has been shown above that the breakup ampli- 
tude is proportional to a half-shell matrix element 
(zo,qo,T,\U\(j)p,q,'El). Substituting, similarly to the cal- 
culation of the elastic scattering amplitude, the finite- 
dimensional operator II s into the expression for the 
breakup amplitude and utilizing the projection rules for 
the free-motion and scattering wave-packets, one gets: 



T* s (p,q) 



T 



Ojo,kj 



>tp£s 
iS(pt) R Ojo,kj 

P*kQj9o ' 

<7o G % > 
q G ij, 
PG A> k , 



(38) 



where p* k — ^J2me* k s and q*- = gfe-i + qj] are momenta 
corresponding to and ij bins respectively and the 
D s k is the momentum width of the A| bin. 

Thus, we just have found that the elastic and breakup 
amplitudes can be calculated directly using the diagonal 
("on-shell") and non-diagonal ("half-shell") matrix ele- 
ments of the same operator^. 

However, some problem still arises here: how to define 
correctly which of the basis states Z k ^ correspond to the 
on-shell states of the Hi, Due to the discretization of 
the spectrum every WP basis state corresponds to the 
energy E k - — e s k * + £* and thus one does not get the 
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exact coincidence of these energies for different "on-shell" 
three-body WP states with the energy E of the initial 
state. In other words, the energy conservation for three- 
body WP states is fulfilled only approximately within 
the corresponding bin widths. To avoid this difficulty we 
apply some energy averaging procedure to the transition 
matrix elements which is quite natural for the lattice-like 
representation. 



Now, using the energy averaged WP amplitudes one can 
get easily the following formula for the breakup ampli- 
tude in the hyperspherical representation: 



47rm K 4 



e «(K) T E S) cos(? = 



(43) 



D. Breakup differential cross section 



C. The energy averaging procedure for the 
breakup amplitudes. 

Let us to rewrite expression pip for the breakup am- 
plitude via wave functions of pair subsystems normal- 
ized to (5-functions in energy. Then, when calculating 
the breakup amplitudes one will need the transition op- 
erator matrix elements of the form (we omit for the sake 
of brevity all the spin labels): 

u{E,e) = (z ,ME-£* )\U\(f>(e)i> (E-s)), (39) 

where |0(e)) is the NN scattering state with energy e, 
\tpo(E— e)) is the wave function of the subHamiltonian 
describing the free motion of the third nucleon relative to 
the NN subsystem, and E is the total energy in the cm. 
system. In the framework of the fully-discretized repre- 
sentation it is quite natural to make an energy averaging 
for the transition matrix elements u(E, e) over the exci- 
tation energy e, which leads to the following integrals: 



U n (E) 



1 



de(z , ME - e* )\U\^(e)ME - e)>, 

(40) 

here {[e„_i, £„]} is some set of intervals, in general inde- 
pendent of the initial partition of the excitation energy 
e. A„ = e n — e rl -i are the corresponding widths. 

Further, by replacing the exact U operator with its 
wave packet counterpart (|32[) and using the projection 
rules for the scattering and free WP's, one can define a 
new approximation for the breakup amplitudes: 

^n S (E) = T o£,^[nun( £n)£fc , E - S^) - (41) 



■ max(e n _i, e fc -i, E - £,•)] . 



Here the sum runs over all possible indices k and j for 
which the difference in the square brackets is positive. 

Then, for the single-component breakup amplitudes, 
one obtains the following approximate expression with 
the elements T„: 



P„ 
S n = 2 \ £ ' 



§m(£ 



-30> 

Sn)- 

(42) 



After the determination of the single-component 
breakup amplitudes, the total breakup amplitude is de- 
rived using contributions of all three Faddeev compo- 
nents and can be written with the help of the particle 
permutation operator P in the following way: 

J 4S t (p,q;q D ) = (p,q|(l + P)iG C/|z ,qo). (44) 

The differential breakup cross section is related to this 
total breakup amplitudes as follows [l|: 

-^f— = (2 7 r) 4 |^fc s K r t (p,q;q )| 2 . (45) 

Here ki and k2 are momenta of particles to be registered, 
S is the arclength of the kinematical curve and k$ is the 
phase space volume defined as: 



«l«-2 



k^.k,,:. k!)\ 2 ( 2fc 1 -k 1 (k lab -k 2 ' x 2 



(46) 



fci 



where ki a b = §qo is the center of mass momentum in the 
laboratory system. 

In the case of s-wave NN interactions, the total 
breakup amplitude is defined as sum of spin-quartet 
(S = |) and spin-doublet (E = |) single-component 
terms defined for all different Jacobi sets. E. g. the 
S'-wave amplitude of the two-neutron emission can be 
represented as a sum of three terms 

Kfl 2 = (2|Mi 1 | 2 + |M5°| 2 + |A/i 1 | 2 )/3, (47) 

where M Ss are the total amplitudes for the quartet and 
doublet channels. In accordance with Ref. [32|], they are 
expressed through single-component breakup amplitudes 
defined for different sets of Jacobi momenta (p a ,q a ) as 
follows (we assume here that neutrons are particles 1 and 
2 and the proton is the particle 3): 

Ml 1 = Ti\p 1 ,q 1 )-Ti 1 (p 2 ,q 2 ), 

Mh ° = 7s{* [Ti0(pu qi) + Ti ° (p2 ' q2}] ~ 

l[Ti 1 {p 1 ,q 1 ) + Ti\p 2 ,q 2 )}+Ti°(p 3 , q3 )}, 
M^= \{T^( Puqi )-T^( p . 2 , q2 ) + 



2 



(P2,q2)+ 

^( Pl , qi )-T^(p 2 ,q 2 )}. 



(48) 
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Finally the differential cross section of the n — d breakup 
is expressed through the partial total amplitude (PfTI) as: 



d 5 a 



tt 2m 



dk x dk 2 dS 4 3<? 



(49) 



To determine the differential cross section for the two- 
neutron emission one has to calculate the elements T„ 
(1411) for every spin component (Ss) of the total breakup 
amplitude and then substitute them into the explicit for- 
mulas dm and (|471) . 

Thus, it has been demonstrated above that in our WP 
approach one can find quite naturally all breakup am- 
plitudes together with the elastic scattering amplitude. 
This gives a very nice universal and unified calculation 
scheme. Some details of the numerical procedure for solv- 
ing the matrix equation (pJ3"| are discussed in Appendix 
C. 



IV. ILLUSTRATIVE EXAMPLES 

Because we have developed here a novel approach for 
the treatment of the three-body breakup processes we 
would need precise and reliable tests to check our new 
procedure. In all the tests below we employ as a conve- 
nient universal WP basis the free momentum wave pack- 
ets for Jacobi momenta q and p constructed using the 
generalized Tchebyshev grid: 



tan 



2N+ 1 



1,...,N, 



(50) 



where g m is the common scale parameter and the t- 
parameter determines the "sparseness degree" of the bin 
set. A similar grid with the size M and common scale 
p m is introduced for discretization of the momenta p. 



A. nd breakup amplitudes for a separable NN 
potential 

As a first extremely convenient test we have chosen the 
three-body model with a separable AW-potential. This 
model, which can be treated numerically very accurately 
in various kinematical breakup situations seems to be 
very appropriate for such a test (see below). 

We consider here the nd breakup with pairwise sepa- 
rable AW-interactions in the form: 



= \ s \tp s )(tp s \, s = 0, 1 



(51) 



As it is well known, the Faddeev equations for such po- 
tentials can be reduced to one-dimensional integral equa- 
tions in momentum space. Such equations, as demon- 
strated in our previous work (24|, can be solved quite 
accurately using the two-body free WP basis. For con- 
venience of the reader we describe all the details of this 



procedure in Appendix B. Here we will refer to the re- 
sults found in such approach as to the "exact ones" . We 
will compare the latter with the results of the solution of 
the general three-body WP scheme with the Eq. ([55]) for 
the separable potential. 




10 15 

e (MeV) 



FIG. 3: (Color online) Comparison of the energy unaveraged 
(solid curve) and averaged (dashed-dotted curve) breakup 
amplitudes T2 1 for the separable AW potential (|51|) cal- 
culated with the general WP technique at the basis size 
M = N= 100 . 



The Fig. [3J shows the unaveraged (i.e. having a his- 
togram form) and energy averaged breakup amplitudes 
for the Yamaguchi potential obtained from the general 
matrix FE ([55]) in the lattice basis. The energy averaged 
amplitudes have been calculated using the averaging pro- 
cedure described above. It is clearly seen from Fig. [3J 
the method developed makes it possible to find rather 
smooth energy dependence for the breakup amplitudes. 
Using the averaged amplitudes, we have found the single- 
component breakup amplitudes T Es as functions of hy- 
perangle 8 for both quartet and doublet three-body spin 
channels. 

Now, we can compare the approximated WP breakup 
amplitudes (|41[) found within our general formalism for a 
separable model with the exact amplitudes derived from 
directly solving the one-dimensional Faddeev equations. 
In the Figs. |4]|6] such a comparison between approximated 
and "exact" results is presented. The two-body free WP 
bases with size M — 200 and N — 200 have been used in 
the calculation of the approximated amplitudes. 

We observe in the Figs. [4j6] quite a very good agree- 
ment between the lattice-approximations and the "exact" 
amplitudes. The corresponding curves are almost indis- 
tinguishable in the figures. The only differences are seen 
at the hyperangle region 9 ~ 90° for the spin-doublet 
T2 amplitude. It is clear also that the WP amplitudes 
calculated in a finite-dimensional Li basis converge to 
the exact ones with increasing basis size. 
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FIG. 4: (Color online) The breakup amplitude T2 1 in the 
quartet channel calculated for the separable NN potential via 
the general WP technique for basis dimensions M = N — 100 
(dotted curves) and M = N = 200 (dash-dotted curves) in 
comparison with the 'exact' values (full curves), with the 
resolution of this figure the three curves can practically not 
be distinguished. 
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FIG. 5: (Color online) The breakup amplitude T2 1 in the 
spin-doublet channel for a separable NN model. The nota- 
tions are the same as in the Fig. [4] 



B. n-d breakup amplitudes for a local NN potential 

Having tested our novel approach using the simple sep- 
arable model for the NN force one can move to a more 
realistic case of a local TV TV-interaction. For this, we 
have chosen the so-called MT I-III NN central poten- 
tial which was frequently used in the past for the test of 
few-body calculations. So, we can compare our results 
for this model with very accurate benchmark calcula- 
tions [33j . For the present WP calculations we use again 
the three-body lattice basis constructed on a Tchebyshev 
two-dimensional grid. 

The results of such a comparison are presented in 
Fig. 0S] for our single-component hyperspherical ampli- 



FIG. 6: (Color online) The breakup amplitude in the 
spin-doublet channel for the separable N TV-model. The nota- 
tions are the same as in the Fig. [4] 

tudes A^ s . One can observe in the Figs. EHH] a quite 
satisfactory general agreement with the results of the 
benchmark calculations [33| except in the region 8 ~ 90°, 
similarly to the case of the separable NN interaction. 
It should be mentioned, that the from the amplitude A 
Eq. (I43[) has an additional factor inversely proportional 
to the relative momentum p compared to the amplitude 
T discussed in the previous subsection. So, the differ- 
ences to the exact solution of A are more visible at the 
region corresponding to small values of p. 
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FIG. 7: (Color online) The breakup amplitudes A? 1 in the 
spin-quartet channel calculated using the WP technique for 
the AW force MT I-III for basis size M = N = 100 (dashed 
curves) and M = N = 200 (solid curves) in comparison with 
the results of the benchmark calculations (solid circles) [33J . 

Similar difficulties at 9 ~ 90° are also observed in 
other works [29rf3l| in which the breakup calculations 
have been done in the configuration space. In Ref. [3l| it 
has been demonstrated that for the correct calculation of 
the breakup amplitudes in the area of small relative NN 
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momenta one has to employ the explicit integral form for 
the breakup amplitude T(p,q). 




FIG. 8: (Color online) The breakup A? 1 amplitudes in spin- 
doublet channel for the MT I-III AW potential. The notations 
are the same as in Fig. [7] 

In our approach the minor disagreement of the WP 
breakup amplitudes and the exact benchmark results at 
very low relative momenta can be related to some un- 
certainties in the determination of the widths and the 
corresponding values of momenta for the WP scattering 
states of the pair continuum very close to the threshold. 
So the case of very low relative momenta in the lattice 
approach deserves a separate study which is under way. 



compared to the results of Ref . [33| . One configuration 
includes the FSI peak while the second one is related 
to the so called "space-star" breakup kinematics. For 
derivation of the Faddeev cross sections we used an in- 
terpolation of the data in the table presented in Ref. [33[ 
for the single-component amplitudes. 
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FIG. 10: (Color online) The differential cross section for two- 
neutron emission at the kinematical configuration: 9i = 45° , 
82 — 50.64°, 0i2 = 180° at the incident neutron energy 
Slab = 42 MeV found with usage of the WP technique (dot- 
ted curve) and the conventional Faddeev calculations (solid 
curve). ("Final state peak".) 
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FIG. 9: (Color online) The breakup amplitudes A^° in the 
spin-doublet channel for the MT I-III NN force. The nota- 
tions are the same as in Fig. [7| 

Using further the single-component amplitudes we 
have found the total (i.e. with inclusion of all three Fad- 
deev components) breakup amplitudes as well as differ- 
ential cross sections of two-neutron emission for different 
kinematical configurations. In Figs. [TO] and QT] the dif- 
ferential cross sections for two-neutron emission with our 
WP technique are presented for two configurations and 
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FIG. 11: (Color online) The same as in the Fig.[l0]but for the 
kinematical configuration 6\ = 53.6°, 62 = 53.6°, 4>i2 = 120°. 
("Space star".) 



In summary, the agreement between conventional and 
lattice results is generally very good. The curves are 
mostly hardly distinguishable in Figs. 1101111 
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V. SUMMARY 

In the present work we generalized the wave-packet 
method developed by the present authors earlier for the 
discretization of the three-body continuum and used it 
for finding the three-body breakup amplitudes. As far 
as the present authors are aware, the study in the work 
is the first where the Faddeev breakup amplitudes are 
obtained completely in the three-body L2-basis. Thus, 
it would be appropriate to enumerate some important 
distinctive features of our lattice-like approach. 

1. Due to projection of the scattering integral equa- 
tions onto the wave-packet L 2 basis corresponding to the 
three-body channel Hamiltonian Hi, we get an explicit 
analytical representation for the three-body channel re- 
solvent Gi that is used in all further calculations. For 
this we employ the version of the integral Faddeev equa- 
tions with the kernel Pv\G\ instead of the conventional 
form PtiGo. This simplifies drastically the whole calcu- 
lations scheme as compared to the conventional one, be- 
cause first we do not need to know the full off-shell pair 
^-matrix at many different energies and in addition we get 
the matrix kernel in a very convenient finite-dimensional 
form. As an input information for the AA-interaction, 
we use only the results of a single diagonalization (for ev- 
ery spin channel) of the AA Hamiltonian matrix. From 
such a diagonalization we get immediately the whole set 
of pseudo-states (the scattering WP's) and partial phase 
shifts at many energies corresponding to these pseudo- 
states [281 ] . 

2. For the matrix of the transition operator, one gets 
an universal linear matrix equation with finite matrix el- 
ements. The diagonal (on-shell) elements of this solution 
determines the elastic scattering amplitudes while the 
non-diagonal (half-shell) elements determines the single- 
component breakup amplitudes up to some known phase 
factor. 

3. The structure of the kernel for the matrix equa- 
tion obtained is very convenient for numerical realiza- 
tion. Due to the fact that the kernel is a product of a 
diagonal matrix, two block matrices and a very sparse 
matrix, it is possible to greatly reduce requirements for 
the RAM storage size and noticeably decrease the com- 
putation time. 

4. The effect of the particle permutation operator in 
the Faddeev kernel is represented now with the help of 
the universal matrix of basis functions overlapping for 
different Jacobi coordinate sets. It allows to avoid very 
time-consuming numerous re-interpolations of the cur- 
rent solutions (at iterations) from one set of Jacobi coor- 
dinates to another one at every iteration step [l|. 

5. Due to an averaging of the integral kernels over the 
cells in momentum space the very complicated energy 
singularities of the kernel above the breakup threshold 
(e.g. the moving branching points etc.) are smoothed 
and one can solve the few-body scattering equations di- 
rectly at real energies, i.e. without any contour defor- 
mation to the complex-energy plane. This fact also fa- 



cilitates enormously the practical solution of few-body 
equations above the breakup threshold. 

6. The comparison of the results obtained in our ap- 
proach with those for the model for the separable AA po- 
tential and with benchmark breakup calculations (with a 
semi-realistic local AA-potential) has demonstrated that 
the WP method allows to get quite accurate three-body 
breakup amplitudes and cross sections. Still, the region 
of very low relative AA momenta requires some addi- 
tional study. Some inaccuracy of our results in this area 
can be related to two factors: (i) a slow convergence of 
the WP amplitudes in this region, and (ii) some uncer- 
tainty of the WP representation of the two-body contin- 
uum in the region of very low relative momenta. We plan 
to devote a special study for the solution of this problem. 

In summary one can conclude that the total lattice- 
like Li discretization of the three-body continuum allows 
to find an accurate solution for the three-body Faddeev 
equations for breakup amplitudes and simplifies enor- 
mously the calculation together with a noticeable reduc- 
tion of the computational cost. 
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Appendix A: The permutation matrix in the lattice 

basis 

In our approach we employ a lattice basis, i.e. a 
basis built by free WP's in momentum space. The 
two-dimensional (three-body) wave packet in momentum 
space are step-like functions of variables p and q: 

(p,q\pmj) = (p,q\®ij) = —^${p € S><)*% & S>j), 



which are nonzero only at the intervals 1>i = [pi~i,pi] 
and = [qj—i, qj] (di and dj are the widths of corre- 
sponding intervals). Such wave packets are normalized 
to unity with the weight dpdq and form an orthonormal 
basis (it is assumed that the intervals are not overlap- 
ping). 

The matrix element of the permutation operator P be- 
tween plane waves has a simple form for the s-wave: 

(pVy£|P|p,g,8E) = Af, s PV, <?>,<?) = 
Af, s 4S( P ' 2 + \q' 2 (p 2 + \q 2 )) 0(1 - M), (A2) 

where Af s , is a spin-channel coupling matrix. The 6- 
function guaranties energy conservation and x is the co- 
sine of the angle between vectors q and q' which (with 
taking into account the <5-function) can be expressed as 
a function of three momenta, e.g. p, q, q' : 

«*,-«■/'. (A3) 

qq 
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The condition |x| < 1 in (|A2I) restricts the allowed values 
of momenta to a region, where the overlap is nonzero. 

To find the matrix elements of the permutation op- 
erator P over the free WPs (|A1|) . one has to integrate 
the function P°(p',q',p,q) over rectangular cells Dij = 
Dj (g) D'^j, = 55-, (g) fl'j, (where the upper prime at 
the interval symbol denote that it refers to a different set 
of Jacobi coordinates) : 



1 



(A4) 



(% j ,\P\V ij ) = - 

J 'didjdi'df 

f f P Q {p\q\p 1 q)dpdqdp'dq'. 
Jd'., ., J® ti 

Actually, the integral in Eq. (|A4[) is reduced to an area 
of two overlapping rectangular areas 23y and ^['f 

Hyperspherical (polar in the s-wave case) coordinates 
are most convenient to calculate such overlaps. Let us 
introduce the reduced (rescaled) momentum variable q: 



q = vmy, 



(A5) 



2ti/3 r 



71/2 ■ 



S 7C/3 f 



7C/6 



R(Q) 



7C/6 



71/3 



7t/2 



a 



2tc/3 



FIG. 12: The integration region in the plane (a, a') which is 
the intersection of the area of allowable values of a, a' — the 
large rectangle S determined by the inequalities (|A10|) — and 
the rectangle R(Q) which boundaries depend on the value of 



then, the energy conservation takes the "homogeneous" 
form p 2 + q 2 = p' 2 + q 2 . The hyperspherical coordinates 



Q, a are introduced as usually: 

q = Qsina, p — Q cos a, Q 2 = p 2 



(A6) 



In these hyperspherical coordinates the integral in 
Eq. (|A4|) takes the following form: 



S(p' 2 + - (p 2 + ^q 2 ))ti(l - \x\)dpdqdp'dq' 

= (4/3) J S(Q 2 - Q' 2 )d{\ - \x\) QdQdaQ'dQ'da' 

= l/3n(%,D' .,), 



(A7) 



where we define the overlapping square: 

n(2)y,2)^,) = J - \x\) d(Q 2 )dada' . (A8) 

Thus we get that the permutation matrix element is di- 
rectly interrelated to this square: 



4 n( %,g^,) 



(A9) 



The condition |x| < 1 can be expressed through the hy- 
perangular variables a, a' as follows: 



7T , 7T 7T 

a \ < a < a 

'3 1 2 1 6 1 



(A10) 



So, the overlap region S(a, a') determined by the condi- 
tion \x\ < 1 is a rectangle in the plane (a, a') restricted 
by four straight lines (see Fig. 



Therefore, the integral in Eq. (IA7|) can be evaluated as 
the external (numerical) integral over Q 2 in the range be- 
tween Q^jn and Q„ lax from the area of intersection of the 
rectangle S and the rectangle R(a min , a max , a' min , a' max ) 
whose vertices depend on Q (see Fig. [P2")l : 



11= / d(Q 2 ) / / dada'. (All) 

The integration limits over Q 2 are equal: 

Qlun = max(pti + <%-i,P')-i + (Qj-i) 2 ), (A12) 
QLx = min( Pl 2 + q 2 ,p' 2 + (q'j) 2 ). (A13) 

If Q^ lax < Q^in then the cells do not overlap and the 
matrix element is equal to 0. 

The coordinates of vertices of the rectangle R(Q) are 
computed directly (see the Fig. [T3] for further explana- 
tions): 

( I ■ 3i— 1 Pi \ 

Oimin\Q) = max(arcsm — — — , arccos — ); 

a max (Q) = min(arcsin ^, arccos (A14) 

a min (Q) — max(arcsm , arccos -4); 

q' p' ■ 

Oi' max (Q) — min(arcsin arccos ^ ). (A15) 

Here, if pi > Q then arccos ^ should be replaced by 0, 
and if <ji > Q then arcsin^- should be replaced by 7r/2; 
the same rule should be applied also to primed values. 
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FIG. 13: On the definition of the integration limits (|A12|) 
- (|A15I) in the variables Q,a,a . The cells Dij and D'ij in 
the plane (p, q) and their polar coordinates (a m i n , a max ) and 



( Q min, a m a x) are shown. 



The area of intersection S n R(Q) in the plane a, a' 
is evaluated analytically by the formulas of elementary 
geometry. 



Appendix B: Wave- packet solution for the 
three-body scattering problem with a separable 
potential 

We consider here a system of three nucleons with equal 
masses m interacting by a separable NN force in s-wave 
spin-singlet (s = 0) and spin-triplet (s = 1) states corre- 
spondingly: 

v, = \.\<p.){<p.\, 8 = 0,1. (Bl) 

The two-body i-matrices are also separable: 

t s {E) = |p,)T, |, (B2) 

where r are known functions 

Tr 1 (E) = K 1 -(<Ps\9o(E)\Vs), (B3) 

and go(E) is the free two-particle resolvent. Function 
t i for triplet state has the pole at the deuteron binding 
energy Sq. The corresponding bound state wavefunction 
\zq) is defined as follows: 



l*o) = ^R(e* )g {e* )\<fi), 



(B4) 



where R(£q) is the residue of r t (E) at the pole. 

We use here the two-parameter Yamaguchi potentials 
with the form factors: 



In this case t s (E) and R(sq) take the form: 

7rm 1 



= A 



R(e* ) 



Ps ((3 S - iyME) 2 

Pl{b 2 +P 2 fp h 



Pb = yj-mel 



(B6) 



The potential parameters /3 S and X s are taken from 
Ref. O. 

The Faddeev equation for the elastic transition 
operator U in this case reduces to the system of 
one-dimensional integral equations of the Lippmann- 
Schwinger type [27] for the elastic scattering amplitudes 
corresponding to the total spin S = ^, 1 and orbital mo- 
mentum L (in the case of the s-wave pair interactions the 
£ and L are conserved separately): 

F* L ( q , q'; E) = Z%{q, q'; E) + £ J (q") 2 dq" X 



ZfL(q',q";E)T sl [E- 



3(g") : 
4m 



-^F^ s ,(q",q';E), (B7) 



where 

F^(q,q';E) = (q, <p„ XL\g (E)U(E)g (E)\q' , <p s ,, EL). 

(B8) 

The kernels ZfJf(q,q';E) in Eq. (|B7| are defined as 
follows: 

Z3W;£) = 

A>f^ t (.) |™^f (B9) 
y_l E — q z /m — q' /m — /m 

where a; = cos(qq'), Pj, are the Legendre polynomials, 
and A^ s , is a spin-channel coupling matrix. In case of 
quartet scattering one has the single equation with s = 1 

3 

and Aj^ = — 1, while there are two coupled equations 

2. 

with s = 0, 1 in the case of doublet scattering and Ag = 
A 2 — 1 A 2 — A 2 — — 2 

iv ll — 2' Jl 01 — ■''■10 ~ 2 

After solving Eq. (IB7|) . the partial wave elastic on-shell 
amplitude can be defined from the relation 



^iol 



^« ii(eo)^u i (9o,go;^), 



(BIO) 



where q = \/ §ra(L — gjj). 

Now we turn to the determination of the breakup am- 
plitude. Substituting the explicit formulas for the t- 
matrix in (|B2I) and the deuteron wave function (|B4I) into 
([26| , one can express the partial "Faddeev" breakup am- 
plitudes via the elastic amplitudes Fltf". 



ip s (p) = (p 2 +P 2 s )- 1 . 



(B5) 



Tf L (p,q) = ^R{^)ip s (p)r s (E-Sq 2 /Am)F]i L (q,qo;E). 

(Bll) 

Let's now proceed with the lattice version for elas- 
tic amplitude. We introduce the free WP basis ([5]) and 
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project Eq. (|B7I) to this basis. Finally, we find the matrix 
equation: 

F = Z + ZtF, (B12) 

where the letters with double lines denote matrices of cor- 
responding operators in the WP subspace for the given 
values of total spin X and orbital momentum L (below 
we shall omit X and L for brevity). More definitely, 

Zf!, = -±= [ dqdq'Z?L(q, q'; E) 

\JdAr 

= "T= I dqdq'F^ (q,q'; E) 

7-; = -^ /_ dqq 2 T s (E-3q 2 /4m)(q,q';E). 
\Jdi Js> * 

The elastic on-shell amplitude in WP representation is 
defined from the diagonal "on-shell" matrix elements of 
the lattice transition matrix X: 

A cl (E)~ 2 ^q / bi£ ^° , E-s* eV j0 (B13) 

Similarly, the packet approximation for the breakup am- 
plitude is determined by off-diagonal matrix elements of 
F: 

, eS , (B14 , 

where T s (q*) — t s (E — 3g* 2 /4m) and q* is the midpoint 
of bin !%. 

The WP representation for the breakup amplitude 
A(9) which determines the asymptotics of the breakup 
wave function in hyperspherical coordinates has the form: 

(B15) 

Appendix C: Features of the numerical procedure 

Here we will discuss some details of the numerical pro- 
cedure for solving the matrix equation (|33p . The main 
difficulty is its large dimension. Quite satisfactory re- 
sults can be obtained with a basis size M ~ N ~ 200. 
It is means that in the simplest one-channel (quartet) 
case one gets a kernel matrix with dimension M x N ~ 



40000 x 40000 which occupies — 6.4 GB (at single pre- 
cision) of RAM or external memory of the computer. In 
the two-channel (doublet) case the required amount of 
memory increases by a factor 4. 

However, the matrix of the kernel K for equation ([33]) 
can be written as the product of four matrices which have 
the specific structure: 

K = PViGi =OP°ViGi, (CI) 

where 

Vi = O t Vi. 

The matrix of the channel resolvent Gi is diagonal and 
its elements are defined by simple explicit formulas. The 
matrix of the potential Vi has a block-type structure 
(fM)) : in fact, it is the direct product of the (M x M)- 
matrix of the two-body interaction and the unit ( N x N)- 
matrix. The rotation matrix O has a similar form and 
the actual dimension (M x M). The free permutation 
matrix P is very sparse due to the energy conservation 
condition. As a rule, only about 1% of its elements are 
distinguished from zero, and the sparsity increases when 
the basis dimension increases. 

So, if to summarize all these details one finds that the 
kernel in the matrix equation (|33|) is the product of four 
matrices: a diagonal one, a very sparse one and two block 
matrices with actual dimension (M x M). If, instead of 
storing the entire matrix K, we store its factors only (for 
the sparse matrix P° we store the nonzero elements only), 
then we can save a huge portion of physical memory: 
for the above example we shall need only ~ 128 MB 
instead of initial 6.4 GB. Such an enormous reduction of 
the required memory allows us to perform calculations 
without using an external memory, which, in its turn, 
reduces the calculation time by approximately one order. 

This possibility to avoid storing a very large amount 
of data is related to the specific procedure used by us for 
the solution of the equation (|3"3")) . As a matter of fact, 
to find the elastic and breakup amplitudes one needs not 
all but only on-shell matrix elements of the transition 
operator. So, each of these elements can be found with- 
out complete solving the matrix equation (|33p but by 
means of a simple iteration procedure with subsequent 
summing the iterations via the Pade-approximant tech- 
nique. If we do not store the entire matrix of kernel K, 
to perform each iteration we need only three additional 
matrix- vector multiplications with the matrices of a spe- 
cial form. 

All these features of the procedure lead to an extremely 
economic calculation scheme which can be realized with 
a usual moderate PC. 
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